Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Focal white matter changes in spasmodic dysphonia: a combined diffusion tensor imaging and neuropathological study.

Spasmodic dysphonia is a neurological disorder characterized by involuntary spasms in the laryngeal muscles during speech production. Although the clinical symptoms are well characterized, the pathophysiology of this voice disorder is unknown. We describe here, for the first time to our knowledge, disorder-specific brain abnormalities in these patients as determined by a combined approach of diffusion tensor imaging (DTI) and postmortem histopathology. We used DTI to identify brain changes and to target those brain regions for neuropathological examination. DTI showed right-sided decrease of fractional anisotropy in the genu of the internal capsule and bilateral increase of overall water diffusivity in the white matter along the corticobulbar/corticospinal tract in 20 spasmodic dysphonia patients compared to 20 healthy subjects. In addition, water diffusivity was bilaterally increased in the lentiform nucleus, ventral thalamus and cerebellar white and grey matter in the patients. These brain changes were substantiated with focal histopathological abnormalities presented as a loss of axonal density and myelin content in the right genu of the internal capsule and clusters of mineral depositions, containing calcium, phosphorus and iron, in the parenchyma and vessel walls of the posterior limb of the internal capsule, putamen, globus pallidus and cerebellum in the postmortem brain tissue from one patient compared to three controls. The specificity of these brain abnormalities is confirmed by their localization, limited only to the corticobulbar/corticospinal tract and its main input/output structures. We also found positive correlation between the diffusivity changes and clinical symptoms of spasmodic dysphonia (r = 0.509, P = 0.037). These brain abnormalities may alter the central control of voluntary voice production and, therefore, may underlie the pathophysiology of this disorder.

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Quantitative MRI of spinal cord injury in a rat model: Correlative studies

Serial quantitative and correlative studies of experimental spinal cord injury (SCI) in rats were conducted using three-dimensional magnetic resonance imaging (MRI). Correlative measures included morphological histopathology, neurobehavioral measures of functional deficit, and biochemical assays for N-acetyl-aspartate (NAA), lactate, pyruvate, and ATP. A spinal cord injury device was characterized and provided a reproducible injury severity. Injuries were moderate and consistent to within $\pm$20% (standard deviation). For MRI, a three-dimensional implementation of the single spin-echo FATE (Fast optimum angle, short TE) pulse sequence was used for rapid acquisition, with a 128 x 128 x 32 (x,y,z) matrix size and a 0.21 x 0.21 x 1.5 mm resolution. These serial studies revealed a bimodal characteristic in the evolution in MRI pathology with time. Early and late phases of SCI pathology were clearly visualized in $T\sb2$-weighted MRI, and these corresponded to specific histopathological changes in the spinal cord. Centralized hypointense MRI regions correlated with evidence of hemorrhagic and necrotic tissue, while surrounding hyperintense regions represented edema or myelomalacia. Unexpectedly, $T\sb2$-weighted MRI pathology contrast at 24 hours after injury appeared to subside before peaking at 72 hours after injury. This change is likely attributable to ongoing secondary injury processes, which may alter local $T\sb2$ values or reduce the natural anisotropy of the spinal cord. MRI, functional, and histological measures all indicated that 72 hours after injury was the temporal maximum for quantitative measures of spinal cord pathology. Thereafter, significant improvement was seen only in neurobehavioral scores. Significant correlations were found between quantitated MRI pathology and histopathology. Also, NAA and lactate levels correlated with behavioral measures of the level of function deficit. Asymmetric (rostral/caudal) changes in NAA and lactate due to injury indicate that rostral and caudal segments from the injury site are affected differently by the injury. These studies indicate that volumetric quantitation of MRI pathology from $T\sb2$-weighted images may play an important role in early prediction of neurologic deficit and spinal cord pathology. The loss of $T\sb2$ contrast at 24 hours suggests MR may be able to detect certain delayed mechanisms of secondary injury which are not resolved by histopathology or other radiological modalities. Furthermore, in vivo proton magnetic resonance spectroscopy (MRS) studies of SCI may provide a valuable addition source of information about changes in regional spinal cord lactate and NAA levels, which are indicative of local metabolic and pathological changes. ^

The role of UV radiation in the development of myeloid malignancy in Rag2 knockout mice

Excessive exposure to the UV radiation present in sunlight can lead to the development of skin cancer in humans. Majority of the UV-induced skin tumors in immune-competent mice are highly antigenic in nature. Additionally, they exhibit a high frequency of mutations in the p53 gene, which arise very early in the course of UV radiation and most of them disappear before the development of skin tumors. ^ Initially, this study was to determine whether UV radiation induces skin tumors much earlier in immune deficient Rag2 knockout mice than in immune-competent mice, and if so, compare their antigenic properties and p53 mutation spectra. However, chronic UV irradiation (10 kJ/m2) induced myeloproliferative disease (MPD) as early as 4 weeks in Rag2 knockout mice instead of skin tumors. Conversely, unirradiated Rag2 knockout mice developed MPD at a low frequency, but the frequency increased with the animal's age. Although the UV-irradiated wild type mice (B6129) developed MPD, its frequency was lower and the occurrence much later than the Rag2 knockout mice. ^ This observation led to our new hypothesis that UV irradiation plays a role in the development of MPD in Rag2 knockout mice. After 4 weeks of UV radiation, both histopathology (myeloid:erythroid ratio, number of blast cells) and flow cytometry (mature myeloid, granulocytes and immature cells) demonstrated an increased number of mice affected with the disease in the UV-irradiated Rag2 knockout group than the other groups. ^ We also investigated the role of cytokines and absence of T and B cells in the development of MPD in the Rag2 knockout mice. Results indicated that IL-3 and IL-3Rα chain expression was upregulated in the spleens of the UV-irradiated Rag2 knockout mice (4 weeks). Reconstitution of the Rag2 knockout mice with T and B cells abrogated the UV-accelerated development of MPD. Both histopathology and flow cytometric analysis (mature myeloid cells, granulocytes) showed a decrease in the number of mice affected with the disease in the UV-irradiated, reconstituted group rather than any other group. In summary, this study provides the first experimental evidence that exposure to UV irradiation can lead to the development of MPD in immune deficient mice. ^

Hypersensitive granulomatous response to mycobacterial glycolipid trehalose 6,6'-dimycolate

Protection against Mycobacterium tuberculosis infection requires an effective cell mediated immune response leading to granuloma formation and organism containment. Trehalose 6,6'-dimycolate (TDM), a glycolipid present on the mycobacterial cell wall, has been implicated as a key component in establishment of the granulomatous response. TDM has potent immunoregulatory and inflammatory properties; the acute response to TDM produces pathology resembling early Mycobacterium tuberculosis infection. We have further developed this model to study TDM-specific cell mediated immune responses that may play a role in the later stages of infection and pathology. Lungs from mice immunized with TDM in the form of a water-oil-water (w/o/w) emulsion demonstrate heightened histological damage, inflammation, lymphocytic infiltration, and vascular endothelial cell damage upon subsequent challenge with TDM. This exacerbated response can be adoptively transferred to naïve mice via transfer of non-adherent lymphocytes from TDM immunized mice. To identify the cell phenotype(s) regulating this response, purified non-adherent cell populations (CD4+ and CD8+ T cells; CD19 + B cells) were isolated from TDM immunized mice, adoptively transferred into naive mice, and subsequently challenged with TDM. Lung histopathology and cytokine production identified CD4+ cells as the critical cell phenotype regulating the TDM-specific hypersensitive response. The role of CD1d in presentation of TDM was examined. CD1d, a molecule known to present lipids to T cells, was identified as critical in development of the hypersensitive response. CD4+ cells were isolated from TDM-immunized CD1d -/- mice and adoptively transferred into naive wild type mice, followed by TDM challenge. These mice were deficient in development of the hypersensitive granulomatous response, signifying the importance of CD1d in the generation of TDM-specific CD4+ cells. The experiments presented in this dissertation provide further evidence for involvement of TDM-specific cell mediated immune response in elicitation of pathological damage during Mycobacterium tuberculosis infection. ^

Lactoferrin: An adjuvant to enhance the efficacy of the BCG vaccine

Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), is a disease with world wide consequences, affecting nearly a third of the world's population. The established vaccine for TB; an attenuated strain of Mycobacterium bovis Calmette Guerin (BCG), has existed virtually unchanged since 1921. Intensive research is focused on developing a TB vaccine that can surpass and improve the existing BCG vaccine. Lactoferrin, an iron binding protein found in mucosal secretions and granules of neutrophils was hypothesized to be an ideal adjuvant to enhance the efficacy of the BCG vaccine. Specifically, Lactoferrin enhanced the ratio of IL-12:IL-10 production from macrophages stimulated with LFS or infected with BCG, indicating the potential to affect T-cell development in vivo. Five different vaccination protocols were investigated for generation of host protective responses against MTB infection using Lactoferrin admixed to the BCG vaccine. Mice immunized and boosted at 2 weeks with BCG/Lactofefrin increased host protection against MTB infection by decreasing organ bacterial load and reducing lung histopathology. The observed postchallenge results paralleled with increasing production of IFN-γ, IL-2, TNF-α, and IL-12 from BCG stimulated splenocytes. In vitro studies examined possible mechanisms of Lactoferrin action on BCG infected macrophages and dendritic cells. Addition of Lactoferrin to BCG infected macrophages and dendritic cells increased stimulation of presensitized CD3+ and CD4+ T-cells. Analysis by fluorescent activated cell sorting (FACS) revealed an increase in surface expression of MHC I and decreased ratio of CD80/86 from BCG infected macrophages cultured with Lactoferrin. In contrast, Lactoferrin decreased surface expression of MHC I, MHC II, CD80, CD86, and CD40, but increased CD 11c, from BCG infected dendritic cells, indicating involvement of adhesion molecules. Overall, these studies indicate that Lactoferrin is a useful and effective adjuvant to improve efficacy of the BCG vaccine by enhancing generation of mycobacterial antigen specific T-cell responses through promotion of antigen presentation and T-cell stimulation.^

Metastatic lung cancer in a new mouse model inheriting p53 and latent K-ras mutations

Lung cancer is a devastating disease with very poor prognosis. The design of better treatments for patients would be greatly aided by mouse models that closely resemble the human disease. The most common type of human lung cancer is adenocarcinoma with frequent metastasis. Unfortunately, current models for this tumor are inadequate due to the absence of metastasis. Based on the molecular findings in human lung cancer and metastatic potential of osteosarcomas in mutant p53 mouse models, I hypothesized that mice with both K-ras and p53 missense mutations might develop metastatic lung adenocarcinomas. Therefore, I incorporated both K-rasLA1 and p53RI72HΔg alleles into mouse lung cells to establish a more faithful model for human lung adenocarcinoma and for translational and mechanistic studies. Mice with both mutations ( K-rasLA1/+ p53R172HΔg/+) developed advanced lung adenocarcinomas with similar histopathology to human tumors. These lung adenocarcinomas were highly aggressive and metastasized to multiple intrathoracic and extrathoracic sites in a pattern similar to that seen in lung cancer patients. This mouse model also showed gender differences in cancer related death and developed pleural mesotheliomas in 23.2% of them. In a preclinical study, the new drug Erlotinib (Tarceva) decreased the number and size of lung lesions in this model. These data demonstrate that this mouse model most closely mimics human metastatic lung adenocarcinoma and provides an invaluable system for translational studies. ^ To screen for important genes for metastasis, gene expression profiles of primary lung adenocarcinomas and metastases were analyzed. Microarray data showed that these two groups were segregated in gene expression and had 79 highly differentially expressed genes (more than 2.5 fold changes and p<0.001). Microarray data of Bub1b, Vimentin and CCAM1 were validated in tumors by quantitative real-time PCR (QPCR). Bub1b , a mitotic checkpoint gene, was overexpressed in metastases and this correlated with more chromosomal abnormalities in metastatic cells. Vimentin, a marker of epithelial-mesenchymal transition (EMT), was also highly expressed in metastases. Interestingly, Twist, a key EMT inducer, was also highly upregulated in metastases by QPCR, and this significantly correlated with the overexpression of Vimentin in the same tumors. These data suggest EMT occurs in lung adenocarcinomas and is a key mechanism for the development of metastasis in K-ras LA1/+ p53R172HΔg/+ mice. Thus, this mouse model provides a unique system to further probe the molecular basis of metastatic lung cancer.^

Prevalence and determinants of colonic polyps in children undergoing colonoscopy: A large hospital-based cross-sectional study

Background. Colorectal polyps are abnormal growths in the wall of the colon including the rectum. The study aims to estimate the prevalence and type of colonic polyps in children undergoing colonoscopic examination at Texas Children's Hospital (TCH) in Houston, Texas during 2000-2007. Also, to examine the factors associated with colonic polyps and the potential determinants of colonic polyps in children undergoing colonoscopy and compare those who had colonic polyps with those who did not on colonoscopy, and determine the significant risk factors of colonic polyps in these children. ^ Methods. We conducted a cross sectional study to analyze data collected at TCH. We obtained demographic, clinical, and histopathology information on consecutive patients who underwent colonoscopy during 2000-2007 from endoscopic records contained in the PEDS-CORI registry (Pediatric Endoscopy Database System-Clinical Outcomes Research Initiative), and abstracted data from the accompanying histopathology reports. ^ Results. We identified 2,693-unique patients, under 18 years of age, who underwent colonoscopy. Approximately 65.5% were white non-Hispanic, and 10.8% African-American. The mean age was 8.7 years and 51.8% were female patients. Polyps were present in 174 patients (6.5%). The most common two histological types were juvenile (60.6%), inflammatory (17.4%). We found that the prevalence of polyps was higher in younger aged children (12.9% in 0-5 years) than in older aged children (4% in 15-17 years), and slightly higher in males than in females (7.9% and 5.4% respectively). For males only, the odds of polyps were statistically significantly higher in Blacks and Hispanics compared to white non Hispanics (OR of 2.2 and 2.1, respectively, and 95% CI of 1.3, 3.9 and 1.3, 3.5 respectively). The indications for colonoscopy were different for children with polyps compared to those without polyps, i.e., 47.0% vs. 19.8% respectively for lower GI bleeding, 2.7% vs. 21.4% respectively for abdominal pain/bloating, and, or 0.9% vs. 9.6% respectively for diarrhea. ^ Conclusion. Colorectal polyps occur in about 1 in 15 children and adolescents undergoing first colonoscopy. The demographic variable of younger age is strongly associated with having polyps irrespective of ethnicity. Lower GI bleeding is strongly related to the presence of colorectal polyps in children and adolescents undergoing colonoscopy.^